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Abstract. Electric polarization phenomena in insulating systems have long been 
described in mean field theory by the (static) Clausius-Mossotti or (dynamic) Lorentz- 
Lorenz polarizabilities. It is here shown, in the strong coupling regime, that a 
thermodynamic phase instability exists in these models. The resulting thermodynamic 
phase diagram coincides with that obtained from Dicke-Preparata model of thermal 
superradiance. 

PACS numbers: 42.50.Fx, 32.10.Dk, 33.15.Kr, 42.65.An 

A property of fundamental importance in the physics and chemistry of condensed 
matter insulators is the dielectric response function 

e = l + 47ix (1) 

In principle, e is largely determined by the polarizability[l-12] a of the constituent 
molecules. Precisely relating e to a is a statistical thermodynamic problem of 
considerable complexity. In the lowest order mean field theory approximation, the 
problem is solved by the Clausius-Mossotti equation [13- 16] 

where v is the volume per molecule. Equivalently, the mean field Clausius-Mossotti 
electric susceptibility is then 

X^(^^^). (3) 



(2) 



1 - {Ana/3v)^ 

Note that the stability of Clausius-Mossotti mean field theory requires the inequality 
for the static polarizability 

Ana < 3v (stable normal phase). (4) 

When the above inequality is violated, the normal phase becomes unstable and a new 
stable phase appears. Our purpose is to show that this new stable phase is nothing but 
the superradiant phase implicit in superradiant [17-26] Dicke-Preparata models. The 
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derivation of the superradiant phase diagram is considerably simplified. To see what 
is involved, one may consider the response functions a, x ^ generalized to finite 
frequency u. 

When the response functions are contemplated at finite frequency, for example in 
the index of refraction n(a;) = ^/£{(JJ) , the Clausius-Mossotti Eqs.(l), (2) and (3) 
become the Lorentz-Lorenz equations [27-30]. These have a surprising feature regarding 
frequency shifts. For example, suppose that a single atom with atomic number Z has a 
strong resonance at frequency a;oo, i-e. with m as the electron mass, 

a{u) ^ (^) , (5) 

Eqs.(3) and (5) then imply 

X(^) « .2 .1 o..... . (6) 



vm J Lo\ — up- 



with a "renormalized" resonant frequency ujq given by 



where a = a{uj — 0) is the static polarizability. The macroscopic Lamb frequency shift 
Ao; = uj^ — ujq is certainly worthy of note. It is due to the coherent collective oscillation 
of a large number of molecular dipole moments. 

Electric dipolar motions in (say) a single atom move at a Bohr frequency uj^ 
determined by the energy difference between two states \i) and |/) connected by a 
dipole matrix element (/|p|«) 7^ 0, i.e. huJoo = (e/ — u)- In a condensed matter 
sample of similar atoms, the dipolar motions in one atom radiate an electromagnetic 
field which can drive a neighboring atomic dipole into oscillation. The motion can 
become collective. In the mean field theory of Clausius-Mossotti and Lorentz-Lorenz, 
the driving collective macroscopic electric field E is related to the dipole moment per 
unit volume P via E = — (47r/3)P. The mean electric field phase locks the oscillations 
of many atoms. The collective motions of the large number of dipoles yields the 
macroscopic Lamb shift frequency renormalization cUoo ^0 in accordance with Eq.(7). 
Mean field E theory thereby implies a new collective Bohr frequency which from a 
quantum mechanical viewpoint arises from more closely spaced macroscopic energy 
levels hwo = {Ef — E^) < hj^ = (ej — e^). 

At the critical temperature Tc for stability, we have from the static polarizability 
stabihty Eq.(4) the condition 

^ 47rQ^(Te) ^ _ ^ (critical point), (8) 

where etc = a(Tc). From the dynamical Eq.(7) viewpoint, the renormalized frequency is 
lowered to zero at the critical temperature \imT-^Tc+{^o{T)} — which again leads to 
Eq.(8). 
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In order to compute a{T) one needs a model Hamiltonian H^oi for a single molecule. 
In the presence of an external electric field F one may then compute the single molecule 
partition function 

$(F, T) = -keT In {Tr e-(^rnoi-P-F)/kBT^j _ 

where p is the electric dipole operator of the molecule. The mean dipole moment is 
determined by 

= -sdT -p-dF, (10) 
and the polarizability of the molecule is given by 

an — lim I --^ I . (11) 

For a "two level atom" corresponding to the Dicke model, one chooses the model 
Hamiltonian 

ff.-p.F=( -t "f ) (12) 



yielding 



$2 = -keT In 2 - keT In { cosh ( \) (13) 



«2(T) = ( ^"l tanh (-^ ) . (14) 



so that 

nJT\ = ( ^ tflnh I 

Employing Eqs.(8) and (14) and the definition of the critical volume per two-state atom 
we find for the critical temperature of the Clausius-Mossotti model for two state atoms 

kBTc2 = (t-J-, —^^^7 TT J • (16) 

The above Eq.(16) is identical to the critical temperature of the Dicke model as derived 
in previous work[31, 32]. 

For atoms in an electronic zero angular momentum (J = 0) state, it is more realistic 
to consider a four state model model with three ( J = 1) degenerate states (Mj = 0, ±1). 
The resulting Hamiltonian four state matrix may then be written as 

/ -A -iiF., -iiFy -iiF,\ 
-fxF^ A 
-I^Fy A 
\-piF, A / 

The free energy implicit in the Hamiltonian matrix of Eq.(17) may be written as 

$ = -A;Brin2 - ksT In \ e'^/^^^ + cosh | V^'+^'l^l' j ^ (ig) 



(17) 
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Figure 1. Shown as a "solid curve" is the critical temperature Tc for the four level 
atom as computed from Eq.(20). The curve divides the (T, v) plane into regions which 
are either superradiant or normal. The "dotted curve" shows the critical temperature 
which follows from the superradiant ISll Isl] Dicke two level atom model. 



SO that 



a{T) 



smh{A/ kBT) 



(19) 



A J exp{-A/kBT) + cosh{A/kBT) ' 

Employing Eqs.(8), (15) and (19) yields the critical temperature Tc for the superradiant 
phase transition found by solving 

sinh(A/A:BT,) ^ / v\ ^ ^^q) 

exp{-A/kBTc) + cosh{ A/ kBTc) \VcJ ' 

The phase diagrams implicit in both Eqs.(16) and (20) are shown in Fig.l. The difference 
between the two level model and the four level model is merely a shift in the phase 
boundary. The existence of the superradiant phase is robust with respect to changes in 
the detailed model structure. 

We have developed a very simple algorithm for finding the regions in the (T, v) 
plane for which the superradiant phase exists. Linear magnetic materials can be either 
paramagnetic or diamagnetic. Linear electrical systems can only occur in nature as 



paraelectrics. Diaelectric systems can be ruled out|Q. It is the absence of diaelectricity 
which accounts for the reliability of the algorithm. The computation of the phase 
diagram requires (i) a reliable computation of the single molecule polarizability, and (ii) 
the dielectric constant of the bulk material. The critical temperature has the form of 
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an implicit equation 



a(Tc) = g{v,Tc) where gcM = 




The Clausius-Mossotti function gcM is vahd in mean field theory. The fiuctuation 
corrections to Clausius-Mossotti mean field theory [34-39], (i.e. the precise g{v,T) 
function) may be computed systematically via perturbation theory. The corrections 
are merely quantitative. Qualitatively, the superradiant phase is shifted but left intact. 
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